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TECHNICAL NOTE 4250 


EFFECT OF THE PROXIMITY OF THE WING FIRST-BENDING 
FREQUENCY AND THE SHORT-PERIOD FREQUENCY ON 
THE AIRPLANE DYNAMIC-RESPONSE FACTOR 
By Carl R. Huss and. James J. Donegan 

SUMMARY 


A study of the effect of the frequency of the lowest wing structural 
mode on the airplane center -of -gravity dynamic-response factor was 
made hy employing simplified transfer functions . It was found that the 
simplified transfer function adequately predicted the maximum value of 
the incremental normal- load-factor response at the airplane center of 
gravity to isosceles triangle pulse elevator inputs. 

The results of the study are presented in the form of preliminary 
design charts which give a comparison between the dynamic-response 
factors of the semirigid case and the airplane longitudinal short-period 
case and between the dynamic-response factors of the semirigid case and 
the steady-state value of the airplane longitudinal short-period 
response. These charts can be used to estimate the first-order effects 
of the addition of a wing-bending degree of freedom on the short-period 
dynamic-response factor and on the maximum dynamic-response factor when 
compared with the steady-state response of the system. The results show 
that a structurally damped frequency greater than six times the short- 
period damped frequency will not affect the dynamic-response factor 
of the semirigid short-period response at the airplane center of gravity 
and that, when the frequencies are equal, the semirigid dynamic-response 
factor may be as much as 1.6 times that of the short period. The results 
also show that the maximum dynamic-response factor can be as much as 2.4 
times the steady-state response of the system, depending upon the ratio 
of the natural frequencies of the structural and short-period modes and 
upon the damping of the two modes. 


INTRODUCTION 


As airplanes have increased in size, speed, and flexibility, anal- 
ysis of the loads, stresses, and deflections associated with the 
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longitudinal short-period mode has become increasingly more complex. This 
complexity results from the need to include not only the aeroelastic 
effects associated with structural deformation but also the dynamic 
effects of structural vibration. Considerable effort is currently being 
expended in the field of dynamic analysis and it has become customary to 
express the dynamic effects of both aeroelasticity and structural vibra- 
tion in terms of a dynamic-response factor which relates the dynamic 
response of the airplane to its steady-state response. The effects of 
flexibility are generally associated with a specific response at the 
center of gravity of the airplane, especially in the preliminary design 
stages; however, these effects at other points on the airframe (such 
as a wing-tip deflection or a strain in a particular structural member) 
are often of interest. 

The present-day use of thin high-aspect-ratio wings on large 
high-speed airplanes has resulted in a lowering of the frequency of the 
wing structural vibratory modes. As a consequence of this reduction in 
stiffness, the frequency of the lowest wing vibratory mode is approaching 
the frequency of the airplane short-period mode. The proximity of the 
frequencies of these two modes has a pronounced effect on the airplane 
dynamic-response factor. Although this effect has been Inawn qualita- 
tively for some time and studies of specif ic_ configurations have been 
made, there has been no simple numerical guide for estimating the 
effects of this design trend. Possibly, this lack is a natural conse- 
quence of the nature of the mathematical transfer functions which 
relate the airplane center-of -gravity response to an incremental change 
in elevator angle. These transfer functions are of a type which is 
usually regarded as being more adaptable to specific studies than to 
gener aliz at ion . 

The purpose of this study was to determine whether the results 
obtained by using the complete transfer functions could also be obtained 
to a high degree of approximation with related but greatly simplified 
transfer functions and whether this simplification was of such a nature 
as to permit generalization of the results. The present paper illustrates 
the nature and validity of the simplification of the transfer function 
used and assesses as to both magnitude and trends the effect of the 
proximity of the frequencies of the lowest wing structural mode and the 
airplane short-period mode on the airplane incremental normal load fac- 
tor at the airplane center of gravity. The results are summarized in 
the form of design charts which, it is believed, will be of value in 
the pre lim inary design stages of an airplane. 



NACA TN ^250 


3 


SYMBOIS 
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t> 


^Zh 

Vi 

a Zh 


‘•eh 


b 

C 0 J C l* 


generalized nondimens ional mass-coupling term between Z 
and h degrees of freedom, a^jpSc. 

generalized nondimens ional mass term of flexible-wing mode 
between elastic wing and h degree of freedom, a^/pSc 

generalized nondimens ional mass-coupling term between 6 
and h degrees of freedom, a 0 h/^ Sc 

generalized mass-coupling term between Z and h degrees 


r b o/2 

of freedom, 2 J -m , w .|f z (y)j - S' 



dy Q , slugs 


generalized mass term of flexible-wing mode between elastic 

pb Q /2 f 2 

wing and h degree of freedom, 2 / Jm' w [f z (y)l 

+ I ' w [f0(y)J 2 | <3y oJ slugs 

generalized mass-coupling term between 0 and h degrees of 

ob /2 

freedom, 2 J ' [-S'/ Z (y) + 1‘^fjjCy) - m' w 2f z (y) + 
S’ w Zf 0 (y)J dy Q , slug-ft 


wing span along elastic axis, ft 

Lmens ional tran 
semirigid case 


. Cy dimensional transfer -function coefficients for 


C'i, C' 2 ) • • • C’^ nondimens ional transfer -function coefficients 

for semirigid case 

Cji force coefficient due to elastic-wing deflection, F^/qS 

C m airplane pitching -moment coefficient about the center of 

gravity, M/qSc 

C N airplane normal-force coefficient at the airplane center of 

gravity, F^ y'qS 
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c 


D = 



wing mean aerodynamic chord, ft 


% 

®P 


f,(y) 


kinetic energy, ft-Xb 
potential energy, ft-lb 
force , lb 

airplane normal force at center, of gravity, positive down- 
ward, lb 

spanwise bending-mode shape along wing elastic axis 


f0(y) 


spanwise twisting-mode shape about wing elastic axis per 
unit tip bending deflection, radians/ft 


g acceleration due to gravity, ft/sec 

H wing-tip deflection, h/c, chords 

h wing-tip deflection of elastic axis due to bending, positive 

downward, ft 

I' section moment of inertia, m’^^, slug-ft 2 /ft 

Kq, K x , K^, Ky, Kg, Ky dimensional rigid transfer ^function coefficients 


K^, K' 


k 

l 


M 

m A> m f’ 


p , K' y , K’o, K' q nondimens ional rigid transfer-function coeffic- 
' y ients 

radius of gyration about Y-axis, chords 
reduced angular frequency, cdc/v 

longitudinal distance from airplane center of gravity to wing 
elastic axis (function of spanwise location), positive 
forward, ft 

pitching moment about center of gravity, ft-lb 
mass, slugs Cr 


m 


w 


section mass, slugs/ft 



incremental normal load factor at airplane center of gravity 
(positive upward), g units 

dynamic pressure, lb/sq ft; also, pitching angular velocity, 
radians /sec 

wing plan-form area, sq ft 

section mass moment about elastic axis, m'^, slug-ft/ft 
Laplace transform variable 

Pit 

period of oscillation, — , sec 

CD 

duration of triangular input, sec 

time, sec 
velocity, fps 

longitudinal displacement, positive forward, ft 
lateral or spanwise displacement, ft 

vertical displacement of airplane center of gravity, positive 
downward, ft 

vertical wing deflection of elastic axis due to wing 
bending, positive downward, ft 

angle of attack positive when wing leading edge is up, radians 
. dimensional quasi-steady transfer-function coefficients 

. . r* nondimens ional quasi-steady transfer-function 

coefficients 

dynamic-response factor at airplane center of gravity, 
( An dyn) max/ /An static 

incremental elevator deflection, positive when trailing edge 
is up, deg 

angle of pitch about airplane center of gravity, positive 
nose up, radians 

sweep angle of elastic wing, deg 
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p nondimens ional airplane mass, rn^ypSc 

i damping parameter, percent of critical damping 

p mass density of air, slugs/cu ft 

0 angle of twist of airfoil in plane perpendicular to elastic 

axis, positive when wing leading edge is up, radians 

a) angular frequency, radians/sec 

Subscripts : 

A airplane 

d damped 

dyn dynamic 

f fuselage; structural 

h flexible -wing degree of freedom 

max maximum 

n natural 

o exposed wing 

sp short period 

sr semirigid 

w wing 

Z vertical degree of freedom 

6 pitching degree of freedom IT 

Dots are used to indicate differentiation with respect to time; for 
example, 6 = d0/dt. The subscripts a, §, £l, h, g, H, and 
6 indicate differentiation with respect to the subscripts; for example. 


'N. 


a 


cLC^i /da • 
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GENERAL CONSIDERATIONS 


In the prelim in ary design stage of an airplane, the designer can, 
with presently available methods, estimate the longitudinal short-period 
dynamic-response characteristics of the center of gravity of a given 
configuration for rigid and quasi-steady airframes. The rigid air- 
frame is defined in this paper as a structure that does not deform or 
vibrate, the quasi-steady airframe as one which can deform but not 
vibrate, and the semirigid airframe as one which can both deform and 
vibrate. The problem that the designer is faced with in this preliminary 
design stage is the effect of the airframe vibratory modes (particularly 
those of the wing since it is usually the most flexible) on the 
quasi-steady airframe longitudinal short-period dynamic response. The 
methods available for calculating these effects are usually rather com- 
plex or require information which would probably not be readily avail- 
able at this stage of the design. The designer needs, therefore, some 
means of estimating these effects which are simple and are based on 
parameters which would be available. 

In this paper such means are presented in the form of preliminary 
design charts which can be used to estimate the effects of the proximity 
of the frequencies of the lowest wing structural mode and the airplane 
short-period mode (quasi-steady case) on the dynamic response at the 
center of gravity of the semirigid airplane. The design charts are 
based on the philosophy that in the preliminary design stage of a par- 
ticular configuration the designer will be able to compute either the 
maximum incremental normal load factor for the quasi-steady case or the 
steady- state value of the incremental normal load factor of the quasi- 
steady case. The charts are restricted to estimating the effects of 
only the lowest vibratory wing bending mode on the incremental normal- 
load-factor response at the center of gravity to elevator control inputs. 
All other structural parts are considered rigid. The charts are further 
restricted to a comparison of dynamic-response factors which are defined 
as 

7 - Ky^/^tatic (1) 

where is the maximum amplitude of the first peak of the 

time history of the incremental normal load factor at the center of 
gravity and ^static is steady-state amplitude of the time history 

of incremental normal load, factor. 
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METHOD OP ANALYSIS 


The procedure followed in this paper for studying the effect of the 
proximity of the frequencies of the lowest structural wing mode and the 
Bhort-period mode on the incremental normal-load dynam i c-response fac- 
tor at the airplane center of gravity was patterned after that of refer- 
ences 1 and 2. Dynamic systems representing the incremental ■ normal- 
load-factor response at the center of gravity to an elevator input and 
defined mathematically by transfer functions were excited by various 
isosceles triangular inputs and the maximum values of the resulting time 
responses were expressed as ratios to the steady-state response factors. 

This procedure was followed for systems having the quasi-steady mode coupled 
with a structural mode ( semirigid case) and for the quasi-steady mode 
alone for a wide range of configurations and frequencies and dampings of 
the two modes. The dynamic-response factors thus obtained for the system 
with two modes were then expressed as ratios to those obtained for the 
system with one mode to determine the effects in question. 

Although triangular inputs were used in this study, it is believed 
that comparable results would be obtained for other shapes of pulse- 
type inputs since the process of expressing the semirigid results as 
ratios to the short-period results tend to eliminate the effects of 
different-shaped inputs. Isosceles triangle inputs were used in this 
paper for the following reasons: they approximate in shape severe 

pilot-imposed inputs; their frequency content could be easily varied by 
changing their duration T^; their frequency content could be made 

sufficient to excite the wing structural mode; and they could be easily 
handled mathematically both by manual calculation and by automatic elec- 
tronic calculation. 

For existing airplanes with high-aspect-ratio, thin, flexible wings, 
the lowest structural frequency is usually associated with wing bending 
and, therefore, wing bending was selected as the lowest structural mode 
for this investigation. The theoretical Bystem chosen for this study 
consisted of three degrees of freedom: freedom in pitch, vertical trans- 

lation, and wing bending. The equations of motion developed by Lagrange's 
method describing this system have been previously established and are 
presented in reference 3* For convenience they are also restated in 
appendix A of this paper. 

The assumptions made in this study included the following: lin- 
earity, no change in airplane forward velocity, small perturbations, 
and rigidity of the fuselage and tail assemblies. These assumptions may 
be summarized by the assumption that the motions of an aircraft with 
flexible wings are described by the equations given in appendix A. It 
was further assumed that the aircraft is statically and dynamically 
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stable longitudinally, that is, that the aircraft short-period mode 
and structural mode are oscillatory and are damped. 

As mentioned earlier in this section, the dynamic systems used in 
this paper were defined mathematically by transfer functions relating 
the incremental normal load factor at the center of gravity to an incre- 
mental elevator input. Some of the terms in the transfer functions 
could be eliminated with small loss in accuracy and the analysis was 
made by using these simplified transfer functions. In order to show 
this relationship, it is first necessary to define the complete transfer 
functions and then demonstrate the simplifications that can be made to 
obtain the simplified but practical transfer f mictions . Hereafter in 
this paper the word "complete" will refer to transfer functions containing 
all the terms and the word "simplified" will refer to the transfer func- 
tion with some of its terms omitted. 


Complete Transfer Functions 

The complete transfer functions relating the incremental normal 
load, factor at the center of gravity to the incremental elevator angle 
input for the semirigid case, the quasi-steady case, and the rigid case 
were obtained from the equations of motion given in appendix A. 

Semirigid case .- The transfer function for the semirigid case 
which defines a system that has both wing quasi-steady deformation and 
wing vibration is, in nondimens ional form. 



4 5 

jS + C'gs + C’^s 


+ cy 5 


+ C l 2 s £ 


2 + C'qs + C’^ 

+ cy + c\ 


(2) 


where the definitions of the C T coefficients are given in appendix B. 
In dimensional form the transfer function may be written as 

C.-S^ + Cs-s^ + C-S 2 + CqS +- C^ 

An^f s ) = 2 6 T o 9 (*) 

A5,y ' if. 3 2 ^ 

e s^ + cy -f- c 2 s* + y + 

where the conversion factors of C 1 to C are given in appendix C. 
The static value of this transfer function is seen to be C^/C^. The 
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characteristic equation may he factored into two quadratic equations 
by Graef fe 1 a method and written as 



CijS^ - 4 * Cgs^ 4- 4- Cgs 4" 


s 2 +-fi| 


sp ' 


Wap 3 + W). 


sp 


s 2 4- 


2{ fW/ + W): 


W 


Quasi-steady case .- As indicated in reference 4 by letting rates 
of wing-tip deflection D 2 h = Dh = 0 in the equations for the semirigid 
case, the transfer function for the quasi-steady case may be formed. 

In this case the wings can deform but do not vibrate. The transfer 
function for the quasi-steady case may be written as 



„ 2 

r^s 4- r^s 4- r 5 

2 

s 4- r^s 4- Tg 


(5) r 


where the r coefficients are defined in appendixes B and C. The static 
value of this f unction is seen to be g. It is interesting to note 

that the static value of the semirigid case is equal to the static value 


of the quasi-steady case C<^/C^ = r^I 


since 


r F 1 

1 0 2 


- C 0 C '4 


and 


P o r, 5 = C o C, 9* 


Rigid case .- By letting D 2 h = Dh = h = 0 in the equation for 
the semirigid case, the transfer function for the rigid case may be 
formed and written as 



s) 


K^S 4- KqB 4- Kr j 


K^S 4- 


K_ 


( 6 ) 
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where the K coefficients are defined in appendixes B and C. Although 
the transfer functions of the rigid case and the quasi- steady case have 
the same form, the transfer coefficients of the two cases differ in that 
the transfer coefficients of the rigid case are modified by the effects 
of wing deformation to give the quasi-steady-case transfer coefficients. 
The static value of this function Is ICt/Ko- 


Simplified Transfer Functions 

In order to determine whether the number of terms in the complete 
transfer functions could be reduced, studies of 15 airplane configura- 
tions were made on an electronic analog computer by us in g the complete 
transfer functions for the semirigid case (eq. (4)) and the quasi-steady 
case (eq. (5))- The configurations used covered a range of wing sweep 
angles from 0° to 60 °, of ratios of wing mass to airplane mass from 0.15 
to 0.50, of airplane center-of -gravity positions from 0 . 25 c to 0.45c, 
and of dynamic pressures from 100 to 800 pounds per square foot. These 
studies indicated that some of the terms in the numerator of the trans- 
fer functions did not contribute appreciably to the max imum value of 
the time history of incremental normal load factor for triangular inputs 
but merely acted as phase shifters and thus were not required for the 
purposes of the present study. Typical results of these studies are 
shown for the sem i rigid case in figures 1 and 2 and for the quasi-steady 
case in figures 3 and- 4. 

Semirigid case .- The contribution of the numerator terms of the 
semirigid transfer function with each numerator coefficient equal to 
unity is shown in figure l(a). In figure l(b) the contribution of each 
of these same numerator terms is shown for typical values of the coef- 
ficients. From plots such as these it is seen that the term makes 

the most important contribution to the maximum value of the incremental 
normal load factor. 

Calculations of incremental normal- load-factor time response to 
isosceles triangle inputs were then made by using only the term 

in the semirigid transfer function. These time histories were compared 
with time histories obtained from the complete transfer function 
(eq. (4)), to determine how well the simplified transfer function 

term only in the numerator) described the maximum value of the time 
response of incremental normal load factor for triangular inputs. 

Typical comparisons are shown in figure 2(a) for the case when the fre- 
quencies of the modes are different and in figure 2(b) for the case 
when the frequencies of the modes are equal. 
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On the basis of such computations it was determined that the complete 
semirigid . transfer function (eq. (4)) could be reduced to 


«, \ __ C 9 _ 

s2 + 2 *f(“n) f s + K)/ 

(7) 

and still adequately describe the maximum value of the time history of 
incremental normal load factor for triangular inputs. The use of the 
word "adequately" in this paper means generally to within about 5 per- 
cent and rarely more than about 10 to 15 percent. 



Quasi-steady case .- A similar procedure was used to dete rmin e the 
contribution of the terms in the numerator of the quasi-steady transfer 
function (eq. ( 5 )) to the maximum value of the time response of the 
incremental normal load factor. In figure 5(a) the contribution of the 
numerator terms of the quasi- steady transfer function with each numerator 
coefficient equal to unity is shown. In figure 5(t>) the contribution 
of each of these same numerator terms is shown for typical values of the 
coefficients. In this case it is seen that is the important term. 

Typical comparison of a time history obtained from the reduced transfer 
function term only in the numerator^ with that obtained from the 

complete transfer function (eq. (5)) Is shown in figure 


From comparisons such as that shown in figure b, it was determined 
that the complete quasi-steady transfer function (eq. ( 5 )) could be 
reduced to 


— (s) = 

^e s 2 + Fjs + r 2 


( 8 ) 


and still adequately describe the maximum value of the time history of 
incremental normal load factor for triangular inputs. 

Semirigid short-period case .- Since the denominator of equation (8) 
does not equal the short-period part of the denominator of equation (j)> 
it was found convenient to define another transfer function. This trans- 
fer function will be called the semirigid short-period case and is 
defined as 
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An 

A5 e 


(s) 


A 


s 2 + 2£ 


bpK) sp e 


s + 



( 9 ) 


where the denominator of equation ( 9 ) is identical to the short-period 
part of the denominator of equation ( 7 ) and the constant A is equal to 



The use of this semirigid short-period transfer function as a 
basis of comparison rather than the quasi-steady case reduced the compu- 
tations to practical proportions. If the quasi-steady case had been, 
used, it would have 'been necessary to estimate a new set of derivatives 
which make up the transfer-function coefficients for each new configura- 
tion and flight condition. However, by defining the semirigid short 
period, it was necessary to choose only the damping and frequency of 
the two modes without regard to the derivatives which determine these 
parameters . 

Actually, the semirigid short-period case is practically equal 
to the quasi-steady case since the damping and frequency of the two cases 
are almost the same for a wide range of configurations and q values. 
(See figs. 5 and 6.) A comparison of the natural frequencies of the 
semirigid short-period case with those of the quasi-steady case for a 
wide range of configurations and q values is shown in figure 5* A 
s imi lar comparison of the damping of the semirigid short-period case 
with that of the quasi-steady case is shown in figure 6. The points 
shown in figures 5 and 6 were computed from the data of reference 3 • 

The data of figure 6 indicate that at the higher values of dynamic 
pressure the damping of the short-period case is greater than that 
of the quasi-steady case. Therefore the maxi mum value of the 
response as well as the maximum dynamic-response factor obtained 
from the short-period case would always be less than that of the quasi- 
steady case. Thus the ratios of maximum dynamic-response factors 
obtained by comparing the semirigid case with the short-period case 
would always be greater than (on the conservative side) or the same as 
those obtained by comparing the semirigid case with the quasi-steady 
case. 


Rigid case . - The complete rigid transfer function (eq. (6)) could 
also be reduced to 
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£S-(s) = 

s 2 + KjS + 


( 10 ) 


in order to define the maximum value of the time history of incremental 
normal load factor for triangular inputs. 


Typical Calculations Illustrating Method 

The method used in this study and the preparation of the desired 
preliminary design charts involved a large "number of calculations and 
plots, typical samples of which are shown in figures 7 arid 8. 

In order to obtain the maximum possible dynamic -response factor ■ 
for the range of the variables, it was first necessary to calculate the 
time response to triangular inputs of varying duration (different fre- 
quency content) for each system (a particular combination of the vari- 
ables). A sample of these calculations is shown in figure 7 for both 
the semirigid short-period and the semirigid cases. Some of these compu- 
tations were carried out on automatic electronic computing equipment, some 
on desk-type computers, and some were carried out by using the tables of 
references 5 and 6 in conjunction with automatic electronic computing 
equipment. The dynamic -response factor defined previously as equation (l) 
was determined for each case by picking the value of the first peak 
of the time histories (see, for example, fig. 7 for T^ = 0.4) and 

dividing it by the static value for the particular system being con- 
sidered. These results were then plotted against the period ratio 
Tf/lgp (ratio of the time base of the input to the natural period of 

the short-period mode) in order to determine the max imum dynamic-response 
factor for each case. A typical plot of this procedure is shown in 
figure 8. The data of figure 8 are for the_same cases as those of fig- 
ure 7, the four points shown in figure 8 having been computed from the 
results shown in figure J. 

From plots such as that shown in figure 8, it was possible to 
ascertain the magnitude and trends of the effects of the proximity of 
the frequencies of the lowest structural mode and the airplane longitu- 
dinal short-period mode on the incremental normal-load dynamic-response 
factor at the airplane center of gravity. The plotting of these calcula- 
tions resulted in the desired preliminary design charts. 
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Range of Variables 

Hie results of this study are believed to be valid over a range of 
variables as follows: dynamic pressure from 100 to 800 pounds per 

square foot, wing sweep angles from 0° to 60 °, ratios of wing mass to 
airplane mass of 0.15 to 0 - 50 ; cent er-of -gravity location from 0.25 
to 0.45 mean aerodynamic chord, ratios of damped wing lowest struc- 
tural frequency to damped airplane longitudinal short-period fre- 
quency from 1 to 15 , and damping of the wing lowest structural mode and 
airplane longitudinal short-period mode from 0 to 95 percent of critical 
damping. 


RESULTS AND DISCUSSION 


The results of this paper are summarized in the form of preliminary 
design charts. As mentioned previously, these design charts were obtained 
from plots such as those of figure 8 covering a complete range of combi- 
nations of the variables. 

The chart given as figure 9 is a plot of the ratios of maximum 

dynamic-response factors f 7 sr ') against the ratio of the 

\ y m sL X / y / max 

structural-mode natural frequency divided by the semirigid short-period 
natural frequency, the damping of the structural mode being held con- 
stant at 2 percent of critical damping. The data of figure 9 were 
reduced to a more compact form by converting the abscissa to a ratio of 
the damped frequency of the structural mode and the damped semirigid 
short-period frequency. This simplification is given as the design 
chart shown in figure 10 . 

The design chart shown in figure 10 indicates that, if the damped 
structural frequency and the damped semirigid short-period frequency 
are equal, the maximum dynamic-response factor of the semirigid case 
will be about 1.6 times the value of the maximum dynamic-response factor 
for the semirigid short-period case. It can also be seen that, when the 
ratio of the damped structural frequency and the damped semirigid short- 
period frequency is greater than about 6 , there is no increase in the 
max imum dynamic-response factor of the semirigid case over the value for 
the semirigid short-period case. 

Increasing the damping of the short-period mode of the semirigid 
case while holding the damping of the structural mode constant also 
results In a decrease in the maximum dynamic-response factor of the 
semirigid case. This effect can be seen from the results shown in 
figure 9 . 
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In order to investigate the effect of structural damping on the 
airplane dynamic-response factor, calculations were made for a semirigid 
short-period mode with a natural frequency of 3.162 radians per second 
and a damping of 38 percent of critical damping coupled to a structural 
mode having variable damping of 0 to 95 percent of critical damping and 
a damped frequency equal to the semirigid short-period damped frequency 
(2.926 radians per second) and equal to 5 times the semirigid short- 
period damped frequency. Dynamic-response-factor ratios for these cases 
are plotted against critical damping of the structural mode in fig- 
ure 11. The result shown in figure 11 indicates that, for a given 
value 'of damping of the short-period part of the semirigid case, an 
increase in the damping of the structural mode results in a decrease in 
the maximum dynamic -response factor of the semirigid case. 

Thus, the data of figures 9 > 10, and 11 indicate that, for a given 
frequency of the short-period mode of the semirigid case, an increase in 
the frequency and/or an increase in damping of the structural mode and/or 
an increase in damping of the semirigid short-period mode result in a 
decrease of the maximum, dynamic-response factor of the semirigid case. 
Figures 9 , 10, and 11 would he the ones used hy a designer in order to 
obtain an estimate of the increase in the incremental normal- load 
short-period dynamic-response factor at the airplane center of gravity 
due to the proximity of the frequency of the lowest wing structural mode 
to that of the airplane longitudinal short-period mode. Use of these 
figures presumes, as mentioned earlier, that the designer would be able 
to estimate the maximum longitudinal short-period response and would 
have an estimate of the lowest wing structural frequency and damping of 
the lowest wing structural mode. 

Another design figure which may be useful is one which gives the 
effect of the proximity of the structural natural frequency to the short- 
period natural frequency on the maximum semirigid dynamic response when 
compared with the semirigid short-period static value. This result was 
easily obtained by plotting the semirigid maximum dynamic -response 
factor for each case (obtained from plots such as fig. 8 and noting that, 
as pointed out earlier, the static value of the semirigid and semirigid 
short-period cases are equal) against the ratio of the structural 
natural frequency to the semirigid short-period natural frequency. Such 
a plot is presented as figure 12. 

The designer could use the chart given in figure 12 under the same 
restrictions as were mentioned for the previous charts with one excep- 
tion. This exception is that in using this chart the designer would 
need to know only the airplane longitudinal short-period steady-state 
response rather than the maximum, short-period response. 
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The effect of the input-time base on the dynamic-response-factor 
ratio may also be of interest and can be determined from plots such as 
that of figure 8. In this case, rather than express the maximum values 
of the dynamic-response factor as ratios, the values of the dynamic- 
response factor of the semirigid and short-period frequencies are 
expressed as ratios at specific values of the period ratio T^/Tgp and 

are plotted against the period ratio. Typical plots of this dynamic - 
response-factor ratio are shown in figure 13 for three values of short- 
period damping. 

The base of the input that gives the maximum dynamic-response 
factor is, of course, different for each case, depending on the damping 
of the two modes. It was usually greater than about 0.7 of the natural 
period of the short period for all the cases studied in this paper. 
Examination of plots, such as those shown in figure 13, indicate that, 
when compared for the same triangle base, the highest ratio of dynamic- 
response factor for a frequency ratio of 1.0 will be obtained from 
triangular inputs with a base equal to 0.6 to 0.8 of the natural period 
of the short period. For frequency ratios greater than 1.0, the ratio 
of dynamic-response factors is greatest for triangles with a base equal 
to less than 0.1 of the natural period of the short period. Thus, it 
^ is difficult to pinpoint a specific triangle base as being the one 

giving the most severe results . 

Since airplanes operate at flight conditions (altitude, airspeed, 
center-of -gravity location) which are constantly changing, the frequency 
ratio for a particular configuration will not be constant. Present- 
day large high-speed airplanes with thin, high-aspect-ratio, flexible 
wings are operating in the frequency- ratio range of roughly 4 to 10. 

The conditions for which the frequency ratio will be a minimum depends 
somewhat on the configuration but, in general, operations at low alti- 
tude, high airspeed, and forward center-of -gravity position should 
result in the lowest frequency ratio. This effect can be seen in fig- 
ure 14 where the effect of dynamic pressure and airplane configuration 
on the proximity of the damped frequency of the structural mode to that 
of the short-period mode is given. The data of this figure were con- 
verted from the data of reference 3 • 

The data of figure l4 show that, for airplanes with unswept wings 
or wings with very little sweep, the frequencies of the modes are brought 
into closer proximity by an increase in the dynamic pressure or, for a 
given dynamic pressure, by moving the center of gravity forward. For 
wings with sweep angles greater than about 30 °, these trends of the 
effects of dynamic pressure and center-of -gravity position on the prox- 
imity of the frequencies of the modes are the same. For these cases, 
however, the aeroelastic effects caused by increasing the dynamic 
« pressure usually cause the short-period mode to become statically 
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unstable (indicated in fig. l4 by the frequency ratio going to infinity) 
before the frequencies of the modes can be brought together. Thus, for 
a given configuration the operating conditions will determine the rela- 
tive proximity of the frequencies of the two modes and at what point on 
the abscissas of the design charts the airplane is operating. 

It is well to emphasize that the preliminary design charts given 
are only meant to give first-order effects and to apply only to systems 
which are statically and dynamically stable. Furthermore, since the 
curve given in figure 10 is an envelope of the maximum values of the 
converted data of figure 9 , it will normally give conservative values 
of the ratio of maximum dynamic-response factors. Finally, for a 
particular design problem a detailed analysis including all the variables 
should be made if the "rule-of-thumb" value for the ratio of maximum 
dynamic-response factors given by the chart indicates the possibility 
of a dangerous situation. 


CONCLUDING REMARKS 


The results of this study of the effect of the frequency of the 
first wing bending mode on the airplane dynamic-response factor indicated 
that the maximum center-of-gravity load-factor response to a triangular- 
shaped pulse elevator input could be adequately determined by using a 
simplified transfer function for the semirigid and quaBi-static cases. 

The use of the short-period part of the sem ir igid transfer function as 
a basis of comparison gave results which were either equal to or on 
the conservative side of those that would have been obtained from the 
quasi- steady transfer function. 

As a result of the reduction in the number of terms obtained by 
using the simplified transfer functions, it was possible to construct 
design charts which provide trends and rule-of-thumb estimates of the 
effect of the frequency of the first wing-bending mode on the airplane 
dynamic -response factor. The charts show that the maximum dynamic- 
response factor for the semirigid case will be 1.6 times that of the 
short-period case when the damped frequencies of the structural mode 
and short -period mode are equal. Furthermore, when the frequency ratio 
is greater than about 6, a lightly damped structural mode has little 
or no effect on the dynamic-response-factor ratio at the airplane 
center of gravity, and, as the damping of the structural mode increases, 
the frequency ratio at which the structural mode has negligible effect 
also decreases. Finally, the charts indicate that the semirigid maxi- 
mum dynamic -response factor can be as much as 2.4 times the steady- 
state value of the system, depending on the damping of the structural 
and short-period modes and on the ratio of the natural frequencies of 
the two modes. 
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The dynamic-response factor for a particular configuration will 
vary with the operating conditions (principally with dynamic pressure) 
hut should be a maximum at flight conditions of low altitude, high air- 
speed, and forward eenter-of-gravity position. 

It should be repeated that for a particular design problem a 
detailed analysis should be made if the rule-of-thumb value given by 
the design charts indicates the possibility of a dangerous situation. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 21, 1958- 
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APPENDIX A 
EQUATIONS OF MOTION 


The equations of motion used in this paper will be given here 
briefly for the convenience of the reader. A complete development of 
these equations is shown in references 3 and 4. The equations are 
derived on the basis of the Lagrangian equation: 


&t\dQ / SQ SQ Q 


(Al) 


kinetic energy 
potential energy 

generalized coordinate 
generalized force 

Hie three generalized coordinates used are: 

Z vertical translation 

0 pitching velocity 

h displacement of wing tip due to bending of elastic wing 

For an unswept wing the flexible-wing mode shape consists of 
bending f z (y); and for a swept wing the flexible-wing mode shape 

consists of bending f z (y) combined with twisting per unit bending 
deflection at the wing tip f^(y). The spanwise bending is usually 
assumed to be parabolic and the spanwise twist, linear. 


where 

E k 

Q 
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The nondime ns io nai equations of motion, thus derived are: 
2nD(a - 0) + 2A Zh D 2 H - C^a - Da|| - 

"(a S) - H (°%) - “(*■) - 

2m% 2 I^8 + 2%/h - C^a - Dx| C ^ C^j - 

H (%) - m 'K) - s."* 


2*llA + 2^D(a - 6 ) + + SA^A - Op a - 

Gu 


D0 




- Cp 


AS. 


AS. 


and by definition 


An _ VfDQ _ Do. \ 
AS e gU& e AS e J 


(A4) 


(A5) 


Si mul taneous solution of equations (A2), (A5), (a 4), and (A5) results 
in the semirigid transfer function given as equation (5)* 
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APPENDIX B 

DEFINITION OF TRANSFER-FUNCTION COEFFICIENTS 

The transfer coefficients used for the analysis are defined in 
this appendix . 

Semirigid Case 

The coefficients for the semirigid case are 

C 0 = A hh " A 0h )( 2u " | C N Da ) “ 2 ( A 6h A Zh C m Da + 4pK Y 2 A Zh 2 ) 

/ 

C ’l = cj ' A[lK Y A Zh^ C F a + “ A hh C N a " ^hh 0 ^ + 2 *% 2 9 f bh “ 

V. 

“®Se) • i °»J + (“Wzh + ^Sa)!" 3 * - I %) + 

A Zh C m Da + 2A 0h^ 2 ^ - £ C N Da ) ( 2A Zh + | C F <1 ) + itA 8h( A 0h C N a ~ 

\ 

^Sx) + ^^cSs. + 2Az ^S > 
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°' 2 - ^ - 1 <\> 2 + - Wz + “*<*) * 

(I \% E - 21 %% E ♦ ^%)(^ - i %J - - 

' 5 ‘VHdf 2 * 3 * + 5 °0 + " “sAa ' 

^Sh " 2 ° P IH Cm Daj(' 2,i ' 5 \j + %pWTt a ' *2b(°B a * 

Cs dh)] + ' VSa + 2 °%Sa) ' S) H ( 2|i ' 

5 + i %)} 

Cj = i |[-*H«/o % - <^(3* - I + c^-an - i + 

*^{%\ - \\) * - \\) + [1 s% - 

+ f %)]( 2 ^ - I - (V^B + 2 S*x % ' 
° II aSj( 2Azh + 2 %) + (° F a 0 “DH ' ^V 211 ' V^IE " 

2 ' I %) + (V%I ' V^ H ' 2 'W*e)5 %} 
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°\ ■ C 0 |[l °N a S + °=u( _2tl ' 2 ' %) + (\\ ' 

V«i)(^ + 5 x) + [v^C* 1 - 1 %)] + 1 sv%j 

c ’ 5 ( kA ^\ - 4 j WWS> - 


’6 = iSj [^h C N 5 ( C iV, + ° m q) " ^“feO^Ea + %) + 2A eh C N & (2A Zh - 
C “ EE ) + 2 ( A Zh C m B “ A bh C N 6 )( 2A Zh + I ^gj + 2C m 5 (^h 0 ®^ “ 
2V 2 ) + !«jVo 


DH N 5 


, = JL 
7 " gc 0 


-4^ Cj^A^k + SAg^Cju^ + C F a ) + 2A tta ^ia^Ng “ 


C R G. 
a 


1 ” 8 ) + 2A ® h (°%S " ^a 0 "®) + ( C “nH°% ‘ ° i, dhS)( £A Z>i 

%) + 2 + %) ‘ 5 VMS* + %) + 


2C 


N s 


^^ 2 9f h “ A Zh. G m I)H ~ ^eh 1 
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C 'e = S 0 0 {v^^Sa + %) - *“^S(S* + %) + S S^H ' 

SSk + I S(Sa + %)] + S[SSn - SSh ' 2 S(S» + 

Si)] + ^(SS ' SS) + (SS " c %S)( 2Azh + 2 s)} 
c, 9 ■ ^ (SS ' ssX 24 !** 2 - s) + S(SS " ss) 

Quasi-Steady Case 

The coefficients used in the analysis for the quasi-steady case are 
• as follows: 

r 0 = 2 ^K y 2 (2h - | C^) ^hh^ 2 ~ CF h ^ + A 9h- C ^> x C % + + 

" 2 %J 

r 'l -T 0 [f^S - I %( 2 “ - 5 Sj + 5 Sa(- 2 " - t S)] - 

S) + Mss - ss )- -Is)- %( 2tI - 

2 ° H D a )( 2Azh + I s) - C N h ( 2M % 2 S + A Zi < \) - I C ^ C %K + 
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r '2 * ^{[2 C N a S + - i C \)]( 2 V l£2 - %) + ( C S a S ■ 

V%) (“a + I %) + % - I %) + I < Wh} 

r v- 4p^%\K k2 - <0 + ^(v% - \\) 

r \ - i^{[! SA - 1 sK* + %)] (^ 2 - \) + 

2A Zh( C i» 8 C K E ‘ %%) + f^Zh + 5 C F q )( 0 “5 0 “ H " Cm & C %)} 

r ’ 5 - - w)^* 2 - %) + c ^{\% - 
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£* 0,1 


c 



+ 



K 1 o = JL 
8 8*0 


- 2 °«0 + Sf' 2 * 1 ' 5 °>J + 2 °N 6 (Sco + \) 


K '9 " gKo {- 2 <\) 


:\ 
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APPENDIX C 

CONVERSION FACTORS FOR D3MENS IONALIZING THE NONDIMENS IONAL 
TRANSFER -FUNCTION COEFFICIENTS 


The conversion factors for dimensionalizing the nondimens ional 
transfer functions for the semirigid case are as follows: 




The conversion factors for dimensionalizing the nondimens ional 
transfer functions for the quasi-steady case are as follows: 



r 


5 



2 


r* 


5 
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The conversion factors for dimers ionalizing the nondimens ional 
transfer functions for the rigid case are as follows : 
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C7» Ce, Cg 



Time, t, sec 

(a) Contribution of numerator terms with, coefficients equal to uni ty. 



Tine, t, sec 

(b) Contribution of numerator terms for typical values of the coeffi- 
cients: C g = 565. 25; Cg = 50.576; C ? = 0.604; Cg = -0-504; and 

c 5 = -0.059^- 


Figure 1.- Contribution of individual numerator terms of the semirigid 

a + Crs^ + C Y s 2 + Cos + C Q 

transfer function ^-( s) = —2- -1 2 Z to the 

s^ + CjS? + C 2 s 2 + C ? s + Cj^ 


complete response to a unit-amplitude isosceles triangle input 
with T ± - 1.0, C 1 = 15-4414, C 2 = H 6 . 858 O, C^ « 550.6659, and 

= 554.5269. 
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(a) Modes apart; ( £D d) f /( aJ d) S p = 



(b) Modes together; j{^ 3 ) = 1*0« 


Figure 2.- Comparison of the response to a unit- amplitude isosceles 
triangle input with - 1.0 obtained from using the simplified 

semirigid transfer function ^ ■ ( s ) = 


A6 C 


S 4* CjS^ + CgS^ + C^S + C]^ 


and the complete semirigid transfer function 


£&-(b) = 


C^s + C^s CyS + Cqs + 


4 „ 

s + C^s 


^ + C 0 s 2 + C,s + Ci 



u A 
2 


H .2 


("b ) Contribution of numerator terms for typical values of the coeffi- 
cients: r 5 = 7*093; = 0.054; and = - 0 . 0555 . 


Figure 3*- Contribution of individual numerator terms of the quasi-steady 


transfer function ==-( s ) = 

A5 e 


r 3 s + V + r 

= — 5 to the complete response 

+ r^s + r 2 


to a unit -amplitude isosceles triangle input with = 1 . 0 , 
r x = 3*4328, and r„ = 6.9934. 


Incremental normal load factor 



Figure 4-.- Comparison of the response to a unit-amplitude isosceles triangle input with 
= 1.0 obtained from using the simplified quasi-steady transfer function 


An 

A5_ 


>) = 


+ rtjs + Tg 


and the complete quasi- steady transfer function 
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Figure 6.- Comparison of critical damping of the semirigid short-period 
mode with that of the quasi-steady mode at a center-of-gravity loca- 
tion of 0.25c. 
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Figure 8,- Typical plot of dynamic-response factor 7 against period ratio T^/T^ for the 
semirigid short-period case and semirigid case with a frequency ratio 

of 1.9. (%) * 3-162; & sp « 0.58; (o^) = 6.0; ^ = 0.02. 
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Figure 9 -~ Effect on the ratio of maximum dynamic-response factors of the proximity of the 
structural natural frequency to the semirigid short-period natural frequency, = 0.02. 
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Figure 10.- Effect on. the ratio of maximum dynamic-response factors of the proximity of the 
structural damped frequency to the semirigid short-period damped frequency. - 0.02. 
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Figure 11.- Effect of damping of the structural mode on the ratio of maximum dynamic-response 
factors for two ratios of the damped structural and, semirigid short-period frequencies for 
semirigid short-period damping of O. 58 . *=■ 
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Figure 12.- Effect of the proximity of the structural natural frequency to the semirigid 
short-period natural frequency on the semirigid marlim™ dynamic-response factor. 
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Ratio of dynamic -response factors Y sr /Y sp for constant-period ratios 
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Figure 15 .- Typical plots of the effect of the period ratio on the ratio 

various frequency ratios. 





Itynamic pressure, q, Ib/sq ft 


Figure l4.- Effect of dynamic pressure q on the frequency ratio 
(^d) ^ OT var ^ OUB wing sweep angles, wing-mass-to- 

airplane-mass ratios, and center-of-gravity locations. 
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